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ABSTRACT 

We present observations of the intermediate to massive star-forming region 
105345+3157 using the molecular line tracer CS(2 - 1) with CARMA to reveal the 
properties of the dense gas cores. Seven gas cores are identified in the integrated in- 
tensity map of CS(2 — 1). Among these, core 1 and core 3 have counterparts in the 
A = 2.7 mm continuum data. We suggest that core 1 and core 3 are star- forming 
cores that may already or will very soon harbor young massive protostars. The total 
masses of core 1 estimated from the LTE method and dust emission by assuming a 
gas-to-dust ratio are 5±1 Mq and 18±6 Mq, and that of core 3 are 15±7 Mq and 
11±3 Mq. The spectrum of core 3 shows blue-skewed self-absorption, which suggests 
gas infall - a collapsing core. The observed broad linewidths of the seven gas cores 
indicate non-thermal motions. These non-thermal motions can be interactions with 
nearby outflows or due to the initial turbulence; the former is observed, while the role 
of initial turbulence is less certain. Finally, the virial masses of the gas cores are larger 
than the LTE masses, which for a bound core implies a requirement on the external 
pressure of ^ 10* K cm^'^. The cores have the potential to further form massive stars. 

Key words: infrared: ISM - open clusters and associations: individual (IRAS 
05345-1-3157) - radio continuum: ISM - radio lines: ISM - stars: formation - tech- 
niques: interferometric 



INTRODUCTION 

Despite t he well-k nown fact that most stars form in clus- 
ters (e.g.. lLada fc Lada 2003 ). little is known about the de- 
tailed formation process of massive stars or star clusters. 
In contrast to the well-studied mechanisms for low mass 
star formation, the difficulty for massive star formation is 
mainly due to a massive star's high radiative power. The 
strong radiation hinders infalling material in the envelope 
from further accreting onto the central source. Two ma- 
jor theories have been proposed to solve the puzzle. Ac- 
cretion can still be accomplished through a high accretion 
rate g enerated by the turbulent envir onment (jMcKee fc TanI 



rate g enerated by the turbulent envir onment (|McKee fc lani 
l2003h and also through disks (e.g. iKrumholz et all bOOSl ) 



with outflows that gr eatly reduce the radiative pressure 
IjKrumholz et al.ll2005l V On the other hand, massive stars 
could also form from compe titive accretion or stellar mergers 
iBonnell et all I200l1 . 120041 ') in high stellar density regimes. 
The formation of star clus ters can be due to fra gmentation 
of molecular clouds (e.g., iBate fc Bonnelj|2005l ). although 
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IKrumholz et alj l|2007h concluded that the fragmentation of 
pre-stellar massive cores can be avoided through heating, 
allowing the formation of a massive star from one massive 
core. Only through a detailed investigation of the initial con- 
ditions in a protocluster system can we better address the 
problem of massive star formation. 

Massive protostars have strong interactions with their 
ambient environment. They disperse their natal clouds 
th rough stellar winds and UV radiation (see the review 
bv lZinnecker fc Yorkell2007^ . After such dynamical interac- 
tions, the ambient environment changes enormously and no 
longer preserves the initial conditions from when massive 
protostars were first born. Therefore, in order to examine the 
very early conditions of the onset of massive star formation, 
identifying stella r objects at early evolutionary stages (e.g., 
ISeale et al.1 20091') in prot oclusters appears to be very essen- 
tial. |Beuther_eLalJ (120071 ) suggested that the early evolution 
has four stages starting with High-Mass Starless Cores (HM- 
SCs), which then accrete low and intermediate mass proto- 
stars, leading to High-Mass Protostellar Objects (HMPOs) 
and finally the (still deeply embedded) final star. 

Here we report on the massive star-forming region 
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Figure 1. Contour map of CS J = 2 — 1 integrated intensity with the FWHM beam size shown in the lower right corner. The dashed circle 
is the size of the primary beam of the 10-meter dishes. The synthesized beam size, shown in the right-bottom corner, is 4.3" X 5. 9" with 
P.A.=-57.5°. The noise is 0.25 Jy beam~^ km s~^, and the contours are 5a, Oct, 7a, 8a, 9a, 10a, 12a, Ha, positive and negative. The 
boxes are drawn based on the 3(T detections and are used to define cores. 



around the FIR-bright IRAS source 05345+3157, hereafter 
105345. The IRAS source is associated with an IR clus- 
ter surrounded by m olecular gas at a distance of 1.8 kpc 
(JHennineet al.lll992l ). HCO+(2-l) emission east of 105345 
presents complex and clumpy structures, suggesting that 
the source is on the verge of collapse through the contrac- 
tion of local density peaks (|Molinari et al. || 2002l). A massive 
outflow has been o bserved in C0(2-l) I Zhang et al.l l2005l : 
iFontani et al.l 20091) . a typ ical signature of young stellar 
objects. iKlein et al.l (|2005|) shows a ring of dust emission 
around the cluster 105345 in 850 ^.m continuum data. Two 
compact and massive dust clumps are located in the ring to 
the east of the cluster: 105345 #1 and #2, wi th the mass of 
48 Mp and 38 M,^, respectively . Futhermore. iFontani et al.l 
l|2008l) and IFontani et al.l (|2009|) concluded several objects 
at different evolutionary stages in this region with SMA 
and PdBI observations, including prestellar core candidates, 
young intermediate to high mass class protostars, and an 
early-B ZAMS star. The same papers also revealed two con- 
densations in N2D"'"(2-1) and the nature of the deuterated 
cores have been studied in detail. 

The previous observations showed the possibility of ac- 
tive star formation in 105345. The main goal of this paper 
is to study star formation by investigating the dense gas in 
the 105345 environment. Our strategy uses the CS(2 — 1) 
transition as a probe of gas cores to comprehensively study 
gas properties. We use the D array configuration of the 
Combined Array for Research in Millimeter-wave Astron- 
omy (CARMA) to resolve-out large-scale structures and 
trace dense gas fragments (~ 10^ cm~'^) inside the cloud. 
In this paper we present ~5" angular resolution data of 
CS (2 — 1) and 2.7 mm continuum. In the companion paper 
bv lKlein et al.l (|201lh (hereafter paper II), we examine the 



infrared properties of the fragments. With the two papers, 
we obtain a general picture of star formation in 105345. 



2 OBSERVATIONS 

2.1 CARMA observations 

Observations of the line transition CS J = 2 — ^ 1 
V = 97.981 GHz ) toward the clumps 105345 #1 and #2 



Klein et al.ll2005h were made with CARMA in 2008 June. 
The observations used CARMA D array at A = 2.7 mm. 
The CS(2 — 1) transition was placed in a correlator win- 
dow of 8 MHz width and a velocity resolution of 0.37 km 
s~^. Two 500 MHz bands were set up for continuum obser- 
vations. The phase center was chosen between 105345 #1 
and #2 at RA = 05''37™52.6" and DEC = -F31°59'48.0" 
(J2000) using the local standard of rest velocity Vlsr of 
-18.6 kms"^ (jSchrever et al.lll996l ). 0555-^398 was used as 
the phase calibrator. 3C454.3 and Uranus were used for the 
bandpass calibration and the fiux calibration, respectively. 
The system temperature ranged from 50 K to 300 K with 
200 K (single sideband) being a representative value. The 
baselines ranged from 11 meters to 150 meters. The chan- 
nel maps were generated with natural weighting, producing 
a synthesized beam of 5.7" x 4.2". The primary beam size 
of the CARMA combined array is 85" at 98 GH z. The data 
reduc tion was done using the MIRIAD package (jSault et al.l 
Il995h . The amplitude calibration is estimated to be 10%, 
and the uncertainties discussed afterwards are only statisti- 
cal not systematic. 
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Table 1 . Parameters of the identified cores 



Label 


R.A. 


Dec. 


ao 


bo 


P.A.o 


a 


b 


P.A. 


R 




(J2000.0) 


(J2000.0) 


(arcsec) 


(arcsec) 


(deg) 


(arcsec) 


(arcsec) 


(deg) 


(10-2 pc) 


Core 1 


05:37:53.03 


31:59:34.9 


12.0±1.6 


6.9 ±0.6 


+83.5 ±5.3 


10.9 + 1.8 


4.9 ±0.8 


78.8 


3.2 ±0.4 


Core 2 


05:37:51.67 


31:59:53.0 


8.8 ±2.2 


7.5 ±1.7 


-25.4 ±55.6 


7.2 ±2.4 


5.5 ±2.3 


-3.7 


2.7 ±0.7 


Cores 


05:37:52.19 


32:00:05.4 


12.3 ±3.9 


8.9 ±1.5 


+63.9 ±27.3 


11.4 ±4.2 


7.0 ±1.9 


58.8 


3.9 ±0.9 


Core 4 


05:37:52.58 


32:00:06.5 


11.8 ±4.0 


7.4 ± 1.2 


-83.6 ±14.9 


10.5 ±4.5 


5.9 ±1.5 


-88.3 


3.4 ±0.9 


Core 5 


05:37:55.14 


31:59:57.0 


12.3 ±2.6 


8.1 ±1.0 


+48.2+17.2 


11.5 ±2.8 


5.9 ±1.4 


45.6 


3.6 ±0.6 


Core 6 


05:37:50.39 


32:00:01.0 


11.8 ±7.9 


6.8 ±1.3 


-14.2 ±20.7 


10.7 ±8.7 


4.6 ±1.9 


-9.3 


3.1 ±1.4 


Core? 


05:37:50.64 


32:00:06.3 


11.1 ±0.9 


7.7 ±0.4 


+37.2 ± 7.7 


10.3 ±1.0 


5.1 ±0.6 


36.0 


3.2 ±0.2 



ao, bo and P.A.q are the observed major axis, minor axis and position angle of the cores, a, b and P.A. are 
the deconvolved major axis, minor axis and position angle. R is the equivalent radius, Sqrtab/2. 



3 RESULTS 



3.1 CS Zeroth Moment Map and Spectra 

Figure [1] shows the velocity-integrated intensity map of 
CS(2 — 1); the dashed circle is approximately the primary 
beam size of the 10-m antennas in CARMA (~ 69") at 
A=2.7 mm. We identify seven cores with peak intensities 
above 5a" (1 ct ~ 0.25 Jy beam"^ km s"^) inside the primary 
beam of the integrated CS map . Due to the interferometric 
effect that introduces the negative components in the map, 
it is possible that the weak emission associated with core 2 
is spurious. The boxes shown in the right panel of Figure 
[l] are used for deriving the parameters of the cores (see the 
paragraph below) and are drawn based on the 3a detection. 

For this study, we only focus on the cores inside the 
primary beam size of the 10-m OVRO dishes and neglect 
the emission outside the primary beam, which has a lower 
sensitivity and signal-to-noise ratio. Core 1 is isolated in 
the south. Cores 3 and 4 are spatially close and entangled. 
In Table \T\ the parameters of each core are specified, in- 
cluding the observed sizes (major axis ao, minor axis bo 
and position angle P.A.o), the deconvolved sizes (major axis 
a, minor axis b and position angle P.A.) and the equiva- 
lent radius R=Va6/2. The major and minor axes are de- 
rived by fitting a circular or elliptical Gaussian profile in 
the velocity-integrated map of CS. At the distance of 1.8 
kpc, the sizes of these cores range from 0.05 pc to 0.08 pc 
(see Table [l]), which is the typical size of transition groups 
from Hot Molecul ar Cores (HMCs) t o Ultracompact HII 
regions (UC HIIs) (IBeuther et al.l 120071 ) or of low- mass en- 
velopes (e.g.. lLoonev et al.ll200d ). 

The spectra of the seven cores averaged over the boxes 
in Figure [T] are shown in Figure [21 Most of the spectra show 
a strong single-peaked emission, while the spectra of core 2 
and 3 are double peaked. The spectra are fitted by Gaussian 
profiles as shown in Figure [21 Since the double peaks in the 
spectra of core 2 and core 3 may result from self-absorption 
(Sect. 14. 1|) . we fit the profiles with one Gaussian emission 
and one Gaussian absorption line superimposed. Table[211ists 
the fitting parameters for each core (amplitude, peak posi- 
tion of velocity, FWHM). Note that the spectrum of core 4 
shows a red wing emission and that of clump 6 shows a blue 
wing emission above the fitted line. The line wings may be 
due to the outfiow driven by core 3 (see Sect. 14.2]) . 



3.2 Millimeter Continuum Data 

Figure [31 shows the contours of the CS intensity map over- 
laid on the 2.7 mm continuum data from the same observa- 
tion. There are three 2.7 mm continuum emission peaks in 
the map (MMl, MM2 and MM3). MMl in the south is the 
brightest core (22.7±7.2 mjy). MM2, the faintest (6.0 ±1.4 
mjy), is in the north accompanied by the brighter MM3 
(13.8 ± 3.6 mJy). These sources are all unresolved. 

Comparing the A=2.7 mm observation with CS cores, 
we see that core 1 is associated with MMl although the 
two peaks are shifted by about 3", smaller than the beam 
size. Core 3 coincides with MM3, while MM2 shows no CS 
emission (see below). There are no other 2.7 mm emission 
features associated with the rest of the CS cores. 

MM3 and M Ml correspond to the clumps 105345 #1 
and #2 for which [Klein et al.l l|2005l ) derived masses of 48 
Mq and 38 Mq, respectively. Cores 2, 3 and 4 are all asso- 
ciated with the northern clump 105345 #1. In addition, the 
three c ores CI, C2, and C3 in 225 GHz continuum observa- 
tion bv [Fontani et al.l (|2008l ) coincide with MM3, MM2, and 
MMl, respectively. 

The cores MMl and MM3 have been consistently de- 
tected in the (sub-)millimeter continuum and also in the 
FIR by Spitzer (see paper II for details). Figure [3 shows the 
Spitzer 24 \im image; the bright northern source corresponds 
to our core 3 and the southern source corresponds to core 
1. It is peculiar that at 2.7 mm the northern source is the 
weaker one of the two while at all other wa velengths from 
24 \ira to 1.3 mm (the 225 GHz observation in lFontani et al.[ 
(2008)) the northern source is the dominant source. We spec- 
ulate that there is extended emission tha t cont ributed to 
the single-dish observation by [Klein et al.l ([2005|), which is 
resolved out by the interferometer. Thus, the fiuxes seen 
by the interferometer are smaller than expected from the 
single-dish observations, and one would expect that the ef- 
fect is more pronounced for the larger core MM3. 

MM2 shows no counterpart at 850 /xm or in the FIR, bu t 
it does have a counterpart at 3.6 cm JMolinari et al.[ [20021 ): 
therefore, the continuum emission from core MM2 is likely 
dominated by free- free emission from an Hll-region with only 
a fraction of the emission from warm dust. 



3.3 Estimating the Core Mass 

We present three approaches to roughly estimate the core 
masses. The first method is the LTE mass calculated from 
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Figure 2. Spectra of identified CS cores. Ylsr is -f8.6 km/s. Tfie dashed lines are the best fit Gaussians. 



the integrated flux of the CS(2 — 1) hne transition. The sec- 
ond method is calculating the virial mass from the line width 
of CS spectra. The third method is based on the 2.7 mm 
continuum assuming dust emission and a dust-to-gas ratio. 



molecular line, a correction factor Ct should be applied (e.g.. 
[Goldsmith fc Langej[l999l ): 



Nt, 



CrN and 



3.3.1 LTE Mass 

Following the stan dard rotation temperature - column den- 
sity analysis (e.g.. lTurneilll991f) . CS(2 — 1) spectra can be 
used to derive the column density, and further the mass 
of the dense gas. Because the transition CS(2 — 1) is not 
optically thin, a correction factor needs to be applied to 
the LTE mass derived from the assumption of an optically 
thin line. For an optically thin molecular line, by assum- 
ing all the transitions are excited by a single temperature 
Tex (LTE) and negligible background continuu m emission, 
the column density ca n be calculated from (e.g.. lMiao et al.l 
ll995l : lMehringerlll995[ ) 



a 



N {cm 



=2.04 X 10^ 



6a0b{arcsec^) 

^v^(GHz^)S^J?{debye■2)gkgI 
X [ / I{Jy beam~ )dv{km s~ )] 



where 6a and Of, are the FWHM of the synthesized beam, Q 
is the partition function, Eu is the upper energy level of the 
transition, y is the line rest frequency, S is the line strength, 
and /i^ is the dipole moment, gi^ and gi are the usual K- 
level and reduced nuclear spin degenerac ies. For our ob- 
serva tion, 61^=5.9", 6ii,=4.3", Q=0.86Tcx ( R ohlfs fc Wilsonl 
I2000I ). Eu=7.05 K, i/=97.980968 GHz, S=2, /i=1.96 debye, 
and gi=gk=i. Here we assume an excitation temperature of 
20 K for dense cores for ming massive stars (see review by 
IZinnecker fc Yorkd 120071 ) . In the case of an optically thick 



1 



where r is the optical depth and can be derived from 



T = - ln[l - 



Tmb 



J(rex)-J(2.73)J 

(JRohlfs fc Wilsonl |2000|), where Tmb is the main beam 
brightness temperature and 

hji 1 



J{T) 



^ gh)i/kT , 



Therefore, by assuming uniform densities, the LTE mass can 
be derived from 

M = fimH^{l.l33abD^) 

where ^ = 2.33 is the mean molecular weight, mn is hydro- 
gen mass, D is the distance (1.8 kpc for 105345), X is the 
abun dance ratio of CS to H2 , which is assumed to be 10~® 
here (jRohlfs fc Wilsonl 120001 ). and a and b are the major 
axis and minor axis of a defined core. Furthermore, we cal- 
culated the number density of molecular hydrogen for these 
cores by assuming ellipsoidal cores and a, b and 6 for the 
three axes: uh^ = Mz,TB/2Mif/(47r/3(o/2)(6/2)^). The col- 
umn density and number density of H2 range from 2.7 x 10^^ 
to 11.6 X 10^^ cm-2, and from 0.6 x lO'' to 1.7 x 10** cm"^, 
respectively. The LTE masses of the different cores range 
from ~ 2 Mq to 15 Mq. Table 3 lists the integrated line 
intensity {It = / Idv), the derived optical depth r, the col- 
umn density of CS (Ncs), the number density of H2 and the 
results of LTE mass. 
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Table 2. Fitting Parameters For Cores 



Label 


Amplitude 


Peak Position 


FWHM 




Jy/beam 


km/s 


km/s 


Core 1 


0.704 ±0.002 


-16.592 ±0.004 


2.626 ±0.009 


Core 2 


2.008 ±0.380 


-18.065 ±0.004 


1.987 ±0.053 




-2.040 ±0.379 


-18.071 ± 0.004 


1.480 ± 0.045 


Core 3 


1.749 ±0.032 


-17.659 ±0.005 


3.237 ±0.018 




-1.524 ±0.031 


-17.441 ±0.003 


1.735 ± 0.016 


Core 4 


0.699 ± 0.003 


-19.894 ±0.006 


2.983 ±0.015 


Core 5 


0.337 ±0.002 


-19.172 ±0.012 


3.940 ±0.029 


Core 6 


0.922 ± 0.047 


-17.832 ±0.021 


0.849 ± 0.049 


Core 7 


0.921 ± 0.003 


-17.600 ±0.003 


1.584 ±0.007 



The smallest core mass that our observations are sensi- 
tive to is ~ 1 Mq . This is the result of the same calculation 
as ab ove assuming a flux e quivalent to 3a in integrated flux 
map. iFontani et al.l ((20091) observed their condensations N 
and S in N2D"'" and N2H"'" and expected them to be low- 
mass cores (< 1 Mq) according to the observed deuterium 
fraction. Thus, it is not surprising that we don't detect their 
condensations N and S. 

In the calculation of the LTE masses, a value of 10"^ for 
the CS(2-1) molecular abundance is assumed and the abun- 
dance is also assumed to be a constant. However, due to C- 
bearing molecules frozen onto dust grains at early stages of 
star formation, a central depletion of CS(2-1 ) is usually seen 
in pr estellar cores (e.g. iTafalla et al.l |2002| : IStahler fc YenI 
1201(3 ') , and therefore the assumption of a constant abundance 
may lead to larger uncertainties in the calculation. Also, ac- 
cording to different models or calculations, t he derived val- 
ues f o r CS(2-1) abundanc e are various (e.g. iBelloche et al.l 
|2002| : IPirogov et al■ll2007^ . In most literatures, the abun- 
dance proflle is described by a central depletion and a nearly 
constant value between 10~^ and 10~* follows at outer parts 
of cores. Since the radii where the abundance drops is usu- 
ally very close to the center of prestellar cores, we simplify 
the LTE mass estimation by only adopting the plateau part 
of the molecular abundance. The choice of 10~® sets the up- 
per limit for the calculation and the derived LTE masses 
could vary by about one order of magnitude depending on 
the chosen values for the abundance. 



3.3.2 Vmal Mass 

By assuming the cores are in virial equilibrium, we can de- 
rive the virial mass using the line width of CS. We assume 
a stan dard isothermal density proflle for t he cores from IShul 
l|l977h , p = por~^ (cf. iLoonev et al.ll2003f ). In this case, the 
virial mass can be derived from 



3.3.3 Total Mass Estimation from Continuum Emission 

Here we adopt the standard technique to estimate the total 
mass from the continuum emission, in our case the 2.7 mm 
continuum data for core 1 and 3 (remember that core 2 
could have considerable free- free emission; Sect. [3^ . This 
method assumes that the continuum emission is from dust 
and a dust-to-gas ratio. By assuming optically thin emission 
and a single temperature (isothermal), the total mass (dust 
+ gas) is 



M - 



Bu(Tdust)n 



where F^ is the flux density, D is the distance to the source, 
R is the gas-to-dust ratio, B^{Tdust) is the Planck function 
at dust temperature Tdust, and k^, is the dust opacity. The 
flux density of core 1 is 22.7 ± 7.2 mjy and that of core 
3 is 13.8 ± 3.6 mJy. For the dust opacity, assuming a gas 
density ~ lO'' cm"^, we adopt the extrapolated value of 0.29 
cm'^ g~^ at A =2.7 mm for prestellar core d ust with thin ice 
mantles from lOssenkopf fc Henningj l| 19941 ). The extrapola- 
tion was done with a power law (k oc A~ ) and a value of 
P = 1.77, obtained by fltting the respective K-values in the 
FIR (A > 300nm). Therefore, with the assumption of 20 K 
for the dust temperature and 100 for the gas-to-dust ratio, 
we derive the mass of 18 ± 6 Mq for core 1 and 11 ± 3 Mq 
for core 3. The derived masses depend much on the assumed 
values for the dust temperature, the opacity, and the dust- 
to-gas ratio. The error of the mass estimates given above 
only reflects the uncertainty in the flux measurement. Tak- 
ing the uncertainties of the assumptions into account the 
estimate is good within one order of magnitude. For exam- 
ple, if we would assume a rather high temperature of 47 K 
(Paper II) for core 3 the core mass would only be 7 Mq. If 
we would use the opacity for dust grains without ice mantels 
the mass would also drop to 7 Mq . 

The calculation shows that the total mass from the LTE 
method and continuum emission of core 1 and core 3 are 
(5±1 Mq, 18±6 Mq) and (15±7 Mq, 11±3 Mq), respec- 
tively. Both methods have uncertainties. For example, the 
results depend on the assumptions of physical conditions 
and quantities, such as the excitation and dust temperature, 
wavelength dependence of the dust opacity etc. In CS(2-1) 
emission, core 3 is brighter than core 1, while it is opposite 
in the 2.7 mm continuum emission (see section ITT] and 13. 2p . 
This result causes the larger LTE mass of core 3 and larger 
mass by assuming a dust-to-gas ratio of core 1. Other mid- 
infrared to FIR observations (see paper II) show a brighter 
core 3 than core 1, consistent with our CS(2-1) observation. 
Therefore, in the following discussion, we use the LTE mass 
to compare i t with the virial mass, a method consistent with 
ISaito et al.l l|2006l ). 



M„i, 



3R AV^ 
'G'81n2 



where AV is the FWHM of the observed molecular line (see 
Table 2). The derived virial mass for each core is listed in 
Table 3. For core 2 and 3 with the double-peak spectrum, 
we use the line width of the emission line to derive the virial 
masses. 



4 DISCUSSION 

4.1 A Collapsing Core 

Figure |4] shows two CS(2 — 1) integrated intensity maps in- 
tegrated over velocities of the blue and red components of 
the self-absorption spectrum of core 3. The boxes drawn in 
these two maps are the same box as in Fig [l] Both of the 
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Figure 3. Left panel: CS(2 — 1) emission overlaid on the gray-scale 2.7 mm continuum data. The contours represent the CS integrated 
intensity with the same contour levels with Fig. 1. Right panel: CS(2 — 1) emission overlaid on the gray-scale Spitzer 24 ]im data. The 
contours arc the same as in the left panel. 



integrated intensities from the blue and red peak are in the 
same box. 

Arguably, this result suggests that the double-peak 
structure comes from the same core at the limit of our reso- 
lution (5"), although we cannot eliminate the possibility of 
two cores along the line of sight. However, the Gaussian fits 
in Table [2] (emission + absorption) describe the line feature 
of core 3 considerably well and provide an explanation to 
the double-peaked spectrum as coming from one core and 
showing self-absorption. Moreover, the abso rption is close 
to the systemic velocity IjFontani et al.ll2009l ') strengthening 
the evidence for self-absorption. However, it is important 
to note that the negative components in the map are not 
causing the self-absorption dip as the negative components 
seen in Fig [T] are not relevant in the box of core 3. To be 
more specific, we carefully examined the channel maps and 
there are only few negative contours in the channel where 
the deepest dip occurs (in the velocity of ~ -17 km/s). Also, 
we mapped with short-spacing data removed and only very 
few negative contours remained, and the deep dip is still 
seen. As a result, the negative contours cannot be the main 
cause of the dip. 

The double-peaked spectrum of core 3 shows one of 
the typical characteristics of gas infall motions: an opti- 
cally thick transition (CS(2 — 1) in our case) produces a 
blue-shifted, asymmetrical profile and central dip, with the 
blue-shifted line brighter than t he red-shifted line (e.g., 
IZhou et al.ll 19931 : IChoi et aLlliggsl ). This double-peaked fea- 
ture is due to a temperature gradient and velocity gradient 
between each layer in a collapsing core (e.g., v{r) ~ j--0-5 
for the inside-out model) with a static envelope producing 
the central self-absorption dip. The blue-shifted line comes 
from the back side of the core and the red-shifted line from 
the front of the core. The lines with higher critical den- 
sity and higher excitation temperature near the center will 
be obscured by the nearby lines with lower critical density 
and excitation temperature in th e red peak, r esulting in a 
stronger blue peak than red peak (|Evang[l999|) . In this case 
of 105345, optically thin lines peak in the dip of the self- 
absorbed lines. The absorption peaks approximately at the 



systemic velocity, also verified by optically thin tracers mea- 
sured bv lFontaniet al.l (|2009l '). 

Similar double-peaked profiles have been observed in 
CS(2 - 1) for d i fferent objects (e.g IChoi et al.l Il995l : 
iTafalla et a"Llll998l : iBelloche et al.ll2002f ). However, the sep- 
aration between the blue peak and the red peak in our ob- 
servation is significantly larger (~ 3 km/s) than the separa- 
tion in these studies (< 1 km/s). This broad separation may 
arise from the larger inf all velocity for inte rmediate mass to 
massive protostars (e.g. IXue fc WtJ 120081 ). The absorption 
feature in the envelope may also affect the separation of the 
two peaks. In addition, similar double-peaked profiles ca n 
be produced by a rotating core IjPavlvuchenkov et al.ll2008l ). 
Detailed fitting with radiative transfer models is necessary 
to better understand other dynamics effects (tur bulence, ro- 
tation etc.) in double-peaked profiles (JBelloche e t al. 200^). 

Although core 2 also shows two peaks in the spectrum, 
the nature of core 2 is less understood. The asymmetry in the 
spectrum of the two peaks is not significant, indicating that 
the gas is less likely to be infalling. The dip is close to the 
systemic velocity, suggesting that the double-peaked feature 
may be caused by the self-absorption, while the fidelity of 
core 2 may be affected by the interferometric effect. 



4.2 Dynamics of Cores 

Compared with the line widths for cores in other low-mass 
star forming regions (e.g.. lOnishi et al.ll2002l ). we observe rel- 
atively large line widths from 1.5 km/s to 3.9 km/s for these 
seven cores. The isothermal sound speed for a region with 
a temperature 20 K is ~ 0.27 km/s, much smaller than the 
observed line widths. Therefore, other non-thermal motions 
play a large role in the observed broad line widths. There 
are two possibili ties for this non-thermal line-broadening 
(jSaito et al.ll2006l ) . One possibility is that these cores formed 
from rather turbulent gas and still are turbulent. The other 
explanation is that the lines are broadened by interactions 
with the outflow or stellar wind from central protostars. 

If initial large internal motions played a role in forming 
cores, higher density is necess ary for a core to bi nd the sys- 
tem with gravitational energy. ISaito et al.l (|2006l ) suggested 
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Table 3. Physical Properties of CS Cores 



Label It 

( Jy/beam km/s) 



(xl0"em-2) (xlO'^cm-3) 



^H2 MlTE 

(gcm-2) (Mq) 



(Mq) (Kcm-3) 



Core 1 


1.94 ±0.01 


0.24 


6.12 ±0.03 


1.28±0.47 


0.20 ±0.00 


5.3 ±1.2 


26.8 ±3.4 


2.28x10* 


Core 2 


1.02 ±1.00 


0.21 


3.17±3.11 


0.59±0.53 


0.11 ±0.10 


2.0 ±2.3 


13.2 ±3.4 


0.84x10* 


Cores 


3.30 ±0.13 


0.46 


11.55 ±0.46 


1.69±1.11 


0.39 ±0.02 


14.8 ±6.8 


49.8 ±11.5 


4.97x10* 


Core 4 


2.19 ±0.01 


0.27 


7.01 ± 0.03 


1.22±0.81 


0.23 ±0.00 


7.0 ±3.5 


37.1 ±9.7 


3.24x10* 


Core 5 


1.39 ±0.01 


0.20 


4.29 ±0.03 


0.74±0.40 


0.14 ±0.00 


4.7 ±1.6 


67.9 ±11.4 


3.46x10* 


Core 6 


0.82 ±0.06 


0.36 


2.73 ±0.20 


0.61±0.70 


0.09 ±0.01 


2.2 ±2.0 


2.7 ±1.3 


0.06x10* 


Core 7 


1.53 ±0.01 


0.38 


5.15 ±0.03 


1.04±0.26 


0.17 ±0.00 


4.4 ±0.7 


9.6 ±0.6 


0.56x10* 
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Figure 4. Maps of CS (2 — 1) integrated intensity of the blue and red part of the self-absorbed line of core 3. Left panel: The map of 
the blue emission peak. The integrated velocities range from -21.6 km s~^ to -17.5 km s~^. The contour levels start from 5(T (it = 0.12 
Jy beam~^) and end at 28.3(t with a logarithmic step of log(\/2cr). Right panel: The map of the right emission peak. The integrated 
velocities range from -17.1 km s~^ to -14.9 km s~^. The contour levels also start from 5cr (a=0.18 Jy beam~^) and end at 28.3 cr with 
a logarithmic step of log(v^a'). The negative contours are not shown for simplicity. 



such a possible correlation between line widths and the aver- 
age H2 density: the larger the line width, the higher the av- 
erage density. Here we did not take the overall line widths of 
core 2 and 3 into consideration since their line widths may be 
affected and broadened by infall motions as discussed above. 
We see a weak trend (density vs. line widths) between core 
1, 6 and 7. The line widths and average H2 density of core 
1, 6 and 7 are 2.63 ± 0.01, 0.85 ± 0.05, 1.58 ± 0.01 km/s and 
(1.28 ±0.47) x 10^ (0.61 ±0.70) x lO", (1.04 ±0.26) x lO** 
cm~"^, respectively. However, considering the large uncer- 
tainties in the data, the positive trend is less obvious and 
it is not clear about the role of the initial turbulence to the 
formation of the cores in this region from our data. 

Core 4 and 5 do not follow the relation b e tween the 
averaged density and line widths. iFontani et al.l (|2009l ) de- 
tected an outflow oriented in the west-east direction possi- 
bly driven by the 226 GHz continuum source Cl-b or Cl-a, 
which corresponds to our core 3. The same paper also in- 
dicated that this outflow is probably interacting with the 
southern portion of the highly deuterated condensation N, 
which corresponds to our core 4, and causes the observed 
broad lines. Therefore, the line width of core 4 is more pos- 
sibly influenced by the interaction with the outflow. The 
red-wing emission of core 4 and the blue-wing emission of 
core 6 may result from the red and blue lobes of the outflow. 



Nevertheless, it is less clear about the mechanism for core 5 
to produce the observed spectrum. 

However, it seems that the gravitational energy is not 
enough to bind the systems. The derived virial masses are 
noticeably la rger than the LTE masses except for core 6 
(see TablelSjI. IWang et alj (|2008l ) studied low to intermediate 
mass cores around MWC 1080 with CS(2-1) and obtained 
similar LTE masses (calculated from the same method as 
described in Section 3.3.1) and virial masses. ISaito et al.l 
(2006) also concluded that non-turbulent cores have a sim- 
ilar virial mass to LTE mass but that the virial masses are 
usually larger than the LTE masses for turbulent cores. The 
line width of core 6 is similar to that of non-turbulent cores 
(~ 0.90 km/s average) and the similarity between its virial 
and LTE mass indicates that this core is bound by gravi- 
tational energy. On the other hand, for cores with a larger 
virial mass than LTE mass, external pressure must be ap- 
plied to maintain a bound core. The virial equilibrium with 
an external pressure (by neglecting magnetic fields and ro- 
tation) can be expressed as: 



= 2U + n- AnW'Pe. 



where \}^\Ma 



{a: velocity dispersion) is the kinetic en- 
is the gravitational energy. The re- 



ergy and Q = -'i^jf- 

sults of Pea; is listed in Table [3] Except core 6, which is 



8 Katherine I. Lee, Leslie W. Looney, Randolf Klein and Shiya Wang 



already gravitationally bound, other cores need an external 
pressure Pea;/k ~ 10* K cm~'^ to help maintain the structure 
of the cores. With the choice of higher molecular abundance 
by one order of magnitude, the external pressure needed to 
maintain the cores would decrease to ~ 10^ K cm~'^. Obser- 
vations have shown that such pres sure is observed in se veral 
massive star-forming regions (e.g.. lMcKee fc Tanll2003l . and 
references therein). Therefore, we suggest that these turbu- 
lent cores with larger kinetic energy may still be bound by 
the external pressure. 

4.3 Implication for Massive Star Formation 

We suggest that core 3 and core 1 are young stellar ob- 
jects. They show many indicators of high-mass star forma- 
tion. Since young stellar objects are surrounded by dust and 
gas, they are not able to be revealed by optical observations 
(and shorter wavelength observations) due to dust extinc- 
tion. Consequently, YSOs are often identified by their IR 
emission excess. The detection i n the Spitzer 24 yim band is 
often an indicator of YSOs (e.g.. lCaulet et aLlbOOST ). We ob- 
served an infrared counterpart of core 3 and 1 in the Spitzer 
image (Paper II), implying that both cores contain warm 
dust, suggestive of a YSO candidate. In addition, the CS 
spectrum of core 3 shows the typical signature of infalling 
gas motion — a double-peak feature with the stronger blue- 
shifted emission than the red-shifted emission. The wide sep- 
aration of the two peaks (~ 3km/s) is probably due to the 
higher masses. 

Gas infall is also suggested by fitting the core's SED 
with a radiative-transfer model in the companion paper (Pa- 
per II). The radiative-transfer modeling indicates that cores 

3 and 1 could be an accreting stellar source with the central 
stellar mass of ~ 5 Mq surrounded by an envelope with a 
mass of order of 10 Mq (mass estimates in this paper for 
core 3: 11 Mq (dust mass), 17 Mq (LTE mass) and core 1: 

4 Mq (LTE mass), 18 Mq (dust mass)). 

Although, the virial masses (20 and 50 Mq for core 1 
and 3, respectively) are large compared to the LTE masses 
(4 and 20 Mq), indicating that the cores may be gravita- 
tionally unbound, the cores can still be bound through the 
support of the external pressure. Therefore, core 3 and 1 are 
possibly accreting the envelope to form a star and poten- 
tially becoming a massive star (M« > 8Mq). 

According to the classificat ion of evolutionary s tages 
of individual high- mass stars in iBeuther et al.l ((20071), the 
formation of massive stars start with High-Mass Starless 
Cores (HMSCs) and then form High-Mass Protostellar Ob- 
jects (HMPOs) via harboring or accreting low/intermediate- 
mass protostars. HMPOs are accreting high-mass protostars 
with masses larger than 8 Mq , and consist of Hot Molecu- 
lar Core (HMC), Hypercompact HII regions (HCHIIs, size 
^ 0.01 pc) and Ultracompact HII regions (UCHII regions, 
size ^ 0.1 pc) at their early phases. Given the sizes of core 1 
and core 3 (0.06 pc ~ 0.08 pc) and the fact that the spectra 
of core 1 and 3 still peak in the far-infrared and have not 
approached near-infrared regimes, we suggest that these two 
cores are possibly on their way to become HMPOs. 

Also, our cores are con sistent with the massi ve star for- 
mation model proposed bv iMcKee fc Tad l|2003l ). The the- 
ory states that turbulent cores in massive star-forming re- 
gions with pressure Pea: — 10* ~ 10^ K cm~^ can be grav- 



itationally bound and form stars with high accretion rates 
(10~^ Mq yr^^) in a short time scale (several times the free- 
fall time). Our cores appear to be turbulent and non-thermal 
with line widths ranging from 1.5 km/s to 3.9 km/s. How- 
ever, while the positive correlation between the line widths 
and average H2 density implies that the initial turbulence 
may also infiuence the formation of these cores, such corre- 
lation is not obvious in our data considering the uncertain- 
ties. For core 3 and 4, the broad line widths may be also 
contributed by the outfiow driven by core 3. With the ex- 
ception of core 6, the cores have a larger virial mass than 
LTE mass. With the help of external pressure in the par- 
ent cloud, these turbulent cores can be bound and have the 
potential to become seeds for collapse in the future. 



5 SUMMARY 

Identifying protocluster members and understanding their 
dynamics are the first step to study the initial conditions 
preserved in the parental clouds and further unveil the pro- 
cess of star formation. In this paper, we present observations 
of the intermediate/high mass star forming region IRAS 
053454-3157. The observations have been performed in the 
fine transition CS(2 - 1) with CARMA D array. At this 
line frequency, we also observed the continuum emission at 
A — 2.7 mm. With our observation, our main conclusions 
are as follows. 

(i) There are seven CS cores identified in the CS moment 
map. Both core 1 and core 3 have counterparts in the 2.7 
mm continuum data (MMl and MM3). 

(ii) The LTE mass of all the cores are calculated based 
on the CS(2-1) data. The LTE masses range from ~ 2 Mq 
to ~ 15 Mq for all the cores. The LTE mass for core 1 is 
5±1 Mq and that for core 3 is 15±7 Mq. In addition, dust 
masses are estimated from the continuum data: core 1 has 
the dust mass of 18 ± 6 Mq and core 3 has the dust mass 
of 11 ±3 Mq. 

(iii) Most of the spectra of the seven cores show a single 
peak. Core 3 shows a double-peaked spectrum with the blue 
emission stronger than the red emission, suggesting infall 
motion of gas. 

(iv) Core 1 and 3 are suggested to be intermediate- to 
high-mass protostellar candidates. 

(v) The linewidths for the cores are larger than the ther- 
mal linewidth at 20 K. The broad linewidth of core 3 and 
4 are probably contributed by the outflow driven by core 3 
(|Fontani et al.li2009 l. However, the role of the initial turbu- 
lence to the formation of the cores is not clear by examining 
the correlation between the linewidths and average H2 den- 
sity with the uncertainties. 

(vi) These cores require an external pressure of ~ 10* K 
cm""^ to keep them bound. Such high press ure is common 
among massive-star forming regions (e.g., iMcKee fc TanI 
I2OO3), suggesting that these cores are possible seeds for fu- 
ture star formation. 
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